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We demonstrate the code-division multiplexed (CDM) readout of eight transition-edge sensor microcalorime- 
ters. The energy resolution is 3.0 eV (full width at half-maximum) or better at 5.9 keV, with a best resolution 
of 2.3 eV and a mean of 2.6 eV over the seven modulated detectors. The flux-summing CDM system is 
described and compared with similar time-division multiplexed (TDM) readout. We show that the y^iVpixcis 
multiplexing disadvantage associated with TDM is not present in CDM. This demonstration establishes CDM 
as both a simple route to higher performance in existing TDM microcalorimetric experiments and a long-term 
approach to reaching higher multiplexing factors. 



The transition-edge sensor-*^ (TES) is an established 
calorimetric detector whose applications include x-ray 
astronomy,^ gamma-ray spectroscopy for analysis of 
nuclear materials,^ and probing molecular dynamics 
through x-ray absorption spectroscopy.'''^ These applica- 
tions demand ever-larger arrays of detectors to increase 
photon throughput. Because TES microcalorimeters are 
operated at sub-Kelvin temperatures, the reduction of 
power dissipation and wire count through multiplexing 
is crucial. Ammeters made from superconducting quan- 
tum interference devices (SQUIDs) are widely used to 
read out TESs due to their low noise, low impedance, 
low power dissipation, and high bandwidth. 

A typical multiplexed array consists of multiple inde- 
pendent amplifier channels each reading out N detec- 
tors. The two most mature multiplexing techniques now 
in use for TESs are time-division multiplexing (TDM)^ 
and MHz-band frequency-division multiplexing (EDM).® 
Neither is ideal for microcalorimetry. In an A^-row TDM 
multiplexer, the SQUID noise aliased into the signal band 
grows as y/N , a consequence of inefficient use of the 
readout bandwidth.^ This noise limits multiplexers for 
high-resolution x-ray or gamma-ray microcalorimeters to 
tens of detectors per amplifier channel. EDM avoids the 
•\/]V noise penalty but has its own limitations when op- 
erated in the MHz range, including physically large filter 
components,^ and the degradation of sensor resolution 
by ac biasing.* 

A third multiplexing technique, code-division multi- 
plexing (CDM),^ is being developed at NIST in two dis- 
tinct configurations: CDM through current summation^° 
(/-CDM), and CDM through flux summation" ($- 
CDM). In this letter, we report on an eight-element array 
of TES microcalorimeters read out through $-CDM. The 
energy resolution at 5.9 keV averaged 2.6 eV full- width 
at half-maximum (EWHM) in the seven modulated de- 
tectors. We summarize the $-CDM design and compare 
its noise with that of a similar TDM system. 

We have fabricated $-CDM multiplexers to read out 



arrays of 4, 8 and 16 detectors. Eigure 1 depicts a four- 
detector array and explains its operation. The TES sig- 
nals are encoded through a Walsh basis set^^ defined 
by the polarity of lithographically patterned inductive 
traces. The four orthogonal combinations of the signals 
are read out in sequence. The encoding matrices used in 
the 4- and 8-dctector designs are 
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where matrix columns represent TES detectors, and the 
rows represent readout rows. The encoding matrix or 
Wg gives the coupling polarity between the signal from 
each detector and each readout row. (Viewed in this way, 
a TDM system uses the identity matrix for encoding: 
Wtdm = I-) The Walsh code switches the polarity of 
each TES but the first, eliminating sensitivity in demod- 
ulated data to any amplifier drift or pickup occurring af- 
ter the modulation (e.g., the 60 Hz power-line harmonics 
visible in Eigure 3b but absent from 3c) . In the most de- 
manding applications, the single unswitched input could 
be used without a TES as a "dark SQUID" noise mon- 
itor. Other than the multiplexer chips, all the hard- 
ware required by (&-CDM (SQUID series arrays, wire- 
bonded cryogenic circuit boards, and room temperature 
electronics^'^) is directly interchangeable with TDM. Ex- 
isting TDM systems thus need no modifications to their 
firmware or to the data acquisition software to be "drop- 
in compatible" with $-CDM. The analysis software must 
be enhanced, however, to demodulate the N channels of 
raw data into the detector timcstreams. 
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FIG. 1. (Color online) A four-row implementation of code- 
division multiplexing by flux summation ($-CDM). The TESs 
are dc-biased and thus on at all times. The current signal from 
TES j inductively couples to all four first-stage SQUID ampli- 
fiers (SQl) with coupling polarity defined by column j in the 
modulation matrix W4 (Equation 1). Oppositely oriented in- 
ductors (red/bold) produce a negative coupling polarity. Each 
row of inductors (shaded boxes) is transformer-coupled to one 
SQl. Rows of SQls are operated with a standard TDM proto- 
col (see Ref. 13): the rows are activated sequentially via /ad*:, 
so the signal from one SQl at a time passes to a second-stage 
SQUID (SQ2). The output of SQ2 is routed to a 100-SQUID, 
series-array amplifier and then to room-temperature electron- 
ics. To keep the three-stage SQUID amplifier in its linear 
range, the multiplexer is run as a flux-locked loop (Ref. 13). 
The series array output (SA-out) is digitally sampled; a flux- 
feedback signal FBI is then applied inductively to each SQl 
to maintain SA-out at a constant value. 



An example of Walsh-eneoded and -deeoded data is 
shown in Figure 2. Four photons arrive during a 20 ms 
window on a four-detector $-CDM array. The top panel 
shows the encoded signal recorded by each first-stage 
SQUID as the four detectors each absorb x-rays. The 
photons strike TES 3, 1, 2, then 4; the encoded pulse 
polarities therefore reflect columns of Equation 1 in the 
same order. The bottom panel shows the reconstructed 
signal currents in the individual TESs over the same 
20 ms. 

Because the encoding matrices arc defined by lithog- 
raphy on the multiplexer chip, details of the inductor 
layout and other on-chip sources of cross-talk produce 
unequal couplings between the detectors and SQUIDs. 
The encoding matrices given in Equations 1 and 2 are 
therefore only idealizations; we have measured the true 
encodings to depart from the ideal at the 1 % to 2 % 
level in 4, 8, and 16-detector CDM multiplexers. Mea- 
surements on a partial 32-detcctor multiplexer suggest 
its non-uniformity will be at the same level. A correction 
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FIG. 2. (Color online) Example raw and demodulated data 
from four detectors in a single 20 ms period, (a) The raw, 
encoded outputs, Rfc, from the SQl in four-detector <1>-CDM 
(with vertical offsets for clarity). The SQl outputs correspond 
to rows 1-4 in Equation 1. Manganese fluorescence x-rays 
struck TESs 3, 1, 2, and 4 at 0, 5, 7, and 11ms. (b) The 
same data demodulated by application of W^^ to show the 
per-detector signal currents. The signal-to-noise is too high 
for the noise to be seen in this example. 



computed in offline analysis"'^^ can reduce this imbalance 
to levels below 0.1%. Corrected decoding matrices are 
used for all $-CDM data in this letter. We find this cor- 
rection to be stable over at least several weeks, allowing 
it to be measured once and then applied in real-time data 
analysis. The demodulation requires computation scal- 
ing as A''^ per data sample for A^ multiplexed detectors. 
Our experience shows that the computational burden will 
not prevent scaling the technique up to at least A^ sa 100. 

A further linear arrival-time correction is applied, 
though it makes no significant difference in the present 
observations. At higher A^, where the time between suc- 
cessive samples grows longer, the correction would help 
by reducing the dependence of demodulated pulse shapes 
on the arrival time. The correction will also be impor- 
tant to reduce cross-talk effects in future data sets with 
larger N and with higher photon rates. 

In a TDM channel of A^ detectors, the wide-band 
SQUID amplifier noise level (at fixed sampling rate) 
scales as -s/a" due to aliasing (Figure 3a) . The relation be- 
tween the detector current noise /Namp(TES) and SQUID 
flux noise $Namp(SQi) in TDM is^ 



^Namp(TES) — *Namp(SQl) 



(3) 



For large A, the coupling mutual inductance Afjn of the 
TES signal to the SQl amplifier must be increased to 
compensate for the higher amplifier noise. Mjn is limited 
on the high end by the dynamic range of the SQl when 
tracking the steep leading edge of photon pulses. 

CDM has the advantage of sampling all detectors at 
all times, while TDM samples each only 1/A^ of the time. 
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Mn Kq fluorescence with 8 CDM detectors 
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FIG. 3. (Color online) The scaling of SQUID-amplifier 
noise in TDM and $-CDM. Noise was measured at 85 mK, 
with TESs superconducting to emphasize the amplifier noise 
(rather than TES noise) at high frequencies. The Johnson- 
noise contribution from the TES shunt resistor dominates be- 
low IkHz. The r = L/R time constant of the shunt resis- 
tance and inductance in the TES bias loop causes the John- 
son noise to roll off above 100 Hz. At high frequencies, the 
SQUID-amplifier noise is dominant. All measurements used 
irow = 640 ns and a 2.5 MHz, one-pole RC filter before the 
digitizer, (a) Noise from a single SQUID channel, referred to 
the first-stage SQUID, when read out with one, two, four, or 
eight TDM rows. Dotted lines show the single-row, high-/ 
noise level (0.37 /^$o/-s/Hz) multiplied by successive powers 
of Due to aliasing, TDM amplifier noise grows with the 
number of rows as \/]V (see Ref. 7). (b) Noise in four- and 
eight-channel CDM readout. The signals, which have not 
been demultiplexed via the Walsh matrix, are referred to the 
first-stage SQUID. Lines are seen at the 60 Hz power line fre- 
quency and its harmonics. The dotted lines show the CDM-4, 
high-/ noise level (0.65 /i"I>o/\/Hz) multiplied by 1 and \/2. 
As in TDM, the aliased SQUID noise scales as y/N. (c) De- 
modulated noise, referred to the TES current, in four- and 
eight-channel CDM. Both approach IQpA/y'Hz at high fre- 
quencies. We omit the unswitched channel from the average, 
making the 60 Hz line no longer visible. 



This means that amplifier bandwidth is used much more 
efficiently in CDM. In practice, this works as follows. 
The encoded SQl signals (Figure 3b) suffer the same 
y/nN multiplex disadvantage as in TDM. In decoding, 
however, TES signals average coherently, while the N 
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FIG. 4. (Color online) Mn Ka x-ray fluorescence spectra mea- 
sured separately by eight TES x-ray calorimeters read out 
with <1?-CDM. Spectra are offset vertically for clarity. These 
data have been analyzed with corrected Walsh codes and a 
linear arrival-time correction, and a Gaussian energy resolu- 
tion has been fit, techniques described previously in Ref. 11. 
All detectors have multiplexed energy resolution better than 
3eV except for TES 1* — the only detector subject to low- 
frequency noise pickup in the SQUID amplifier chain. The 
$-CDM resolution matches or exceeds that found with equiv- 
alent TESs read out by TDM. 



samples of amplifier noise average incoherently. The de- 
modulated TES signal and amplifier noise are therefore 
independent of multiplexing factor (Figure 3c). $-CDM 
thus allows Min to remain low at large A'' without in- 
creasing the amplifier noise as referred to TES current. 
Low mutual inductance, in turn, increases the effective 
dynamic range of the SQUID amplifiers and makes the 
system more robust with fast TESs and in the face of 
high pulse rates. 

Figure 4 shows Mn Ka fluorescence spectra measured 
by eight TES x-ray detectors read out with $-CDM. Afl 
detectors (besides the unswitched TES 1) achieved 3.0 eV 
FWHM energy resolution or better at 5.9 keV. The mean 
resolution of 2.6 eV is better than the best previous mul- 
tiplexed TES measurement at this energy.^ Count rates 
in these data are low (approximately 5 Hz per detector), 
but we anticipate operating with much higher rates in 
the near future. 

We view the demonstration of $-CDM presented here 
as important for two reasons. First, $-CDM chips are 
drop-in compatible with existing 32-row TDM systems 
but have higher performance. They offer an immedi- 
ate path to the kilopixel-scale arrays of high-resolution 
TES microcalorimeters desirable in applications like syn- 
chrotron science and the proposed Athena satellite. Sec- 
ond, the $-CDM system provides a bridge to the even- 
tual development of /-CDM, in which the rapid alter- 
nation of SQUID switches replaces transformer windings 
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as the mechanism for encoding TES signals. An /-CDM 
muhiplexer could scale to hundreds of detectors per am- 
plifier channel,® eventually enabling even megapixel-scale 
arrays. 
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